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We report on the measurement of angular coefficients and the forward and backward asymmetry of Drell- 

Yan dielectron pairs from pp collisions at = 1.96 TeV. The angular distributions are studied as a function 
of the transverse momentum of dielectron pair and Aft, is measured using the event counting method. The 
Lam- Tung (Aq — A2 = 0) relation which is only valid for a spin-1 description of the gluon is also tested. 



1. Introduction 



' The general expression pj for angular distribution of the final state electron in the Collins-Soper(CS) frame 
is given by: 
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Here, 9 and (f> are the polar and azimuthal angles of the electron in the CS frame. The angular coefficients, 
Aq to At, are in general functions of the invariant mass Mu, rapidity y, and transverse momentum Pt of the 
dilcpton in the lab frame Q- When integrated over the differential cross section reduces to: 

cx (1 + cos^ 9) + -Aq{1 - 3 cos^ 6) + A4 cos 9 (2) 

d cos 2 

When integrated over cos9, the differential cross section reduces to: 



^ cx 1 -f cos (p + cos 2(j) + fi-j sin (j) + fi^ sin 2(f> (3) 
d(p 

Calculations which are based on perturbative QCD (pQCD) make definite predictions for all of the angular 
coefficients. For pp J*/Z -> e+e^ X, the angular coefficients A5, Aq and Aj are close to zero [l|, and if 
we integrate over positive and negative y, the angular coefficients Ai and A3 are small. We can use Eq. [2] to 
extract Aq and A4 (averaged over y) and Eq. |3]to extract A2 and A3 (averaged over y) under the assumption 
that A^ and are zero (as is theoretically expected). 

The angular distribution of the final state electron in the quark- antiquark (qq) annihilation process qq — > 
7*/Z — > e+e^ (LO) can be written as 

(X [1 + cos^e) + B cos9. (4) 



dcos9 



Here Afi,{Mu) = |i? is the forward-backward asymmetry (which originates from the 'y* /Z interference). 

In Quantum Chromodynamics (QCD) at the order of as (NLO) this occurs either through the annihilation 
process with a gluon (g) in the final state (qq — )■ j* /Z g), or via the Compton process with a quark in the final 
state {qg — 7> 1*/Z q), as shown in figure [TJ For the qq — >■ j*/Z g annihilation process [3-I3]j NLO predicts 
that the angular coefficients Aq and A2 are equal, independent of Parton Distribution Functions (PDFs), or 

y, and are described by A?!' = A3? = ^ a ■ For the qg — ?► 7*/^ q Compton process, Aq and A2 depend 

on PDFs and y. However, at NLO, when averaged over y, Aq and A2 are approximately [8|, |9| described by 

Aq* — Air' fa j^^^^pr- The equality A2 = Aq is known as the Lam- Tung relation At LO, it is valid 
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for both qq and gq processes Fixed order perturbative QCD calculations [l| at NLO, as well as QCD 
resumniation calculations to all orders indicate that violations of the Lam- Tung relation are small. The 
Lam- Tung relation is only valid for vector (spin 1) gluons. It is badly broken for scalar (spin 0) gluons [ll]. 
Therefore, confirmation of the Lam- Tung relation is a fundamental test of the vector gluon nature of QCD and 
is equivalent to a measurement of the spin of the gluon. Here, we measure the angular coefficients, Aq, A2, A3, 
and A4 as a function of Pt and the first test of the Lam- Tung relation at large dilepton mass and high transverse 
momentum. This measurement provides a detailed test of the production mechanism of gauge boson with finite 
Pt- In addition, we measure Aft {Afb{Mu) — §^4) in mass which is sensitive to weak mixing angle, sin^ 9w- 
In the end, we compare the measurement with pythia and various predictions which are QCD resummation 
calculations and QCD calculations in the fixed order perturbation theories. 




Annihilation diagram Compton diagram 



Figure 1: Leading order annihilation {qq — > "/* /Z g) and Compton {qg — > "/* /Z q) diagrams for the production of Z 
bosons with finite transverse momentum. 



2. Data Sample and Event Selection 

The data sample consists of 2.1 fb^^ collected by CDF II Detector at Fermilab [T2i| during 2004-2007 for 
the measurement of the angular coefficients. Charged particle directions and momenta are measured by an 
open-cell drift chamber (COT), a silicon vertex detector (SVX), and an intermediate silicon layer in a 1.4 T 
magnetic field. Projective-tower-geometry calorimeters and outer muon detectors enclose the magnetic tracking 
volume. The coverage of COT tracking in pseudorapidity is |?7| < 1.2. The energies and directions of electrons, 
photons, and jets are measured by two separate calorimeters: central (jryl < 1.1) and plug (1.1 < \rj\ < 3.6). 
Each calorimeter has an electromagnetic compartment with a shower maximum detector followed by a hadronic 
compartment. 

The events come from the inclusive single high px central electron or the dielectron trigger. The electron 
trigger efficiencies as a function of Et are measured using the data. The overall trigger efficiencies are measured 
to be almost 100%. [H Q 

The dielectron data sample consists of three different e+e pair topologies. The 
CC topology has both electrons in the central calorimeter with \r]\ < 1.1. The CP topology has one electron in 
the central calorimeter with \ri\ < 1.1, and another in the plug calorimeter with 1.2 < \rj\ < 2.8, a good fiducial 
region requirement. The PP topology has two electrons in the plug calorimeter with 1.2 < \r]\ < 2.8. Events 
with at least one electron candidate with Et > 25 GeV for CC and PP events, and Et > 20 GeV for CP events, 
are selected. The second electron candidate is required to have Et > 15 GeV for CC, Et > 25 GeV for PP, 
and Et > 20 GeV for CP events. The asymmetry Et cut in CC topology increases the acceptance in high Pt 
region. For all topologies, the pp collision vertex along the proton direction (z-axis) is required to be within 60 
cm of the center of the detector. 

For the CC topology, oppositely charged electron pairs are required. One electron has tight and another 
has loose central electron ID cuts with COT track requirement [lal- For the CP topology, the central electron 
passes the tight central electron ID cuts and the plug electron passes the standard plug electron ID cuts with 
a silicon track requirement [l^. For the PP topology, both electrons pass the standard plug electron ID cuts 
[isj . In addition, the pseudo-rapidity must have the same sign and both electrons must have a silicon track for 
the PP topology. The track requirement in both legs of all topologies minimizes the background contamination 
down to ~ 0.5%, especially 7 + jet background as well as dijet background. After the selection cuts, we find 
about 140,000 events for all (CC-I-CP-I-PP) topology. 

For Afb measurement, we use twice more data, 4.1 fb^"'^, with the inclusive single high pT central electron 
trigger sample only. The event selection is same with the measurement of the angular coefficients except that 
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Et > 25 GeV for both legs of CC topology is required and PP topology is excluded in the measurement. The 
PP topology has a large charge fake rate, so it is not optimal for the Afb measurement. 



3. Simulation Sample 

The effect of the acceptance on the angular distributions is modeled using the default pythia [l^ Monte Carlo 
(MC) generator combined with a GEANT simulation of the CDF detector [T^]. The pythia generator used 
at CDF has additional ad-hoc tuning [l^ in order to accurately represent the Z boson transverse momentum 
distribution measured in data. Further tuning was introduced in order to ensure that the MC correctly describes 
the rapidity, as well as the correlations between rapidity and transverse momentum that are observed in the 
data. To reconstruct the simulated events in the same way as data, the calorimetry energy scale and resolution, 
and event selection efhciencies (electron ID and tracking) are tuned using data. 



4. Backgrounds 

Due to the two track requirement in the selection, the background level is very small in the sample. The 
remaning background is measured using data and MC simulations and subtract to get the pure signal. The 
background from QCD process (mainly, QCD dijet) is measured using data. 7-l-jet process is also possible 
background for Drell-Yan process, but the track requirment on both legs rejects most of 7-l-jet background. The 
QCD background is measured using the isolation extrapolation method and the mass spectrum fitting method 
from data directly, fl^ The electron is very isolated object, so the isolation energy shape is a good variable to 
distinguish the electron from the jet object which has a broad and flat isolation energy shape. In the isolation 
extrapolation method, we fit the isolation distribution for both signal and background contributions and then 
extrapolate the background from the high isolation tail into the signal region. To measure the background in Ft, 
the background shape versus Pt is obtained from the background sample (data) and its overall normalization 
is set using the background rate measured by the isolation extrapolation and mass spectrum fitting method. 

The other background source is the background from Electroweak process. The background from Electroweak 
process is obtained using MC simulations. The Electroweak background processes considered are WW, WZ, 
inclusive tt, inclusive W-|-jets, and Z— >■ rr process. We generate the MC simulations for each process and 
estimate the background for corresponding integrated luminosity. 

The overall QCD background is 0.30% and the overall Electroweak background is 0.19%. The background in 
P-T for QCD and Electroweak process is shown in Figure [51 All backgrounds are subtracted as a function of the 
angular distributions, cos 9 and </) for the angular coefficients measurement and mass for Afh measurement. 




|-2 I 

10 20 30 40 50 60 70 

Z pt in all y 



Figure 2: The Pt distribution of dielectron pairs. The event are selected in the full Z boson mass region, 66 < Mee < 116 
GeV/c^. The black crosses correspond to data, the blue histogram to CDF pythia MC. The red and green histogram 
correspond to QCD and Electroweak background, respectively. 
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5. Extraction of Angular Coefficients in Pt 

The Collins-Soper 2] CM frame lepton decay angles, cos and 0, are defined by: 



COS0 = „J ^^^2 - (5) 



tan^=^^!±^.^ (6) 
A,. . Qt 

where M is the mass of the dielectron pair, (i?2) is the four- momentum of electron (positron), Pt is the 
transverse momentum of the electron-positron pair, £^ corresponds to -^(^^ ± ^'^), A-' = ^-^ — Z^, Qt is a 

transverse unit vector in the direction of Qt, and Rt is a transverse unit vector in the direction of Pa x Q- 

To extract the angular coefficients in Pt , the data and the simulated MC (pythia) arc binned in five Pp bins 
which are (0,10), (10,20), (20,35), (35,55), and (55 and above) using the reconstructed dielectron Pp. For each 
Pt range, the data and the MC simulated events are binned in cosO and (j), respectively. The MC events are 
re- weighted to generate the expected experimental distributions for a range of values of Ao and A/^ for cos^, and 
A2 and A3 for (j) distribution. Eq. [5] and [3] are used to re- weight the event to change the angular distributions in 
the generated level. A maximum log-likelihood comparison of the data to MC in cos 9 and is used to extract 
the best values of the angular coefficients, (^0,^4) from cos6' and (^2,^3) from 0, that describe the data. To 
extract the angular coefficients in QED born level (before QED radiation), we reweight the events using the 
angular coefficients obtained before QED radiation. 



5.1. Systematic Uncertainties 

Systematic uncertainties are determined for the background estimation, the energy scale and resolution, the 
electron identification efficiency, the silicon tracking efficiency, the boson Pp and rapidity modeling, and material 
modeling. The systematic uncertainties are considered in cosO and 4> distribtuion for Pr bins. An investigation 
of all of systematic errors listed above shows that the uncertainties on the extracted angular coefficients are 
dominated by the statistical errors. Figure |3] compares the level of systematic vs. statistical errors for the 
angular coefficients in Pt. For the total error, the statistical and systematic errors are added in quadrature. 
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Figure 3: The systematic vs. statistical uncertainty for the angular coefficients. The red bars correspond to the statistical 
errors and the blue bars to the systematic errors. The systematic errors are very small compared to the statistical errors. 
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5.2. The Extracted Angular Coefficients in Pt 

The angular coefficients, Aq, ^3, and A^, are extracted from tlie data in tfie mass region, 66 < Mee < 116 
GeV/c^ for Px bins shown in Figure H] and [S] The measurement is compared with the various predictions 
which are PYTHIA [l^, PYTHIA+Ijet [U, QCD calculations in the fixed order perturbation theory (dyrad [19| . 
MADGRAPH ^2d\, POWHEG [2l|, and FEWz) [22^, and QCD calculations with resummation (resbos [23*1 and 
VBP 24]). The pythia and vbp predictions are close in Pt and follow the approximation of qq annihilation 

p2 

process, -httttt- Other higher order predictions like dyrad, madgraph, powheg, fewz, and PYTHiA+ljet 
are close and have the higher Aq and Ai values than PYTHIA in high Z Pp region. The resbos prediction 
follows PYTHIA in low Pt region, but cross over to the high order predictions in Pt > 35 GeV. 

Aq and Ai have a strong Pt dependence. The measured Aq and A2 follow pythia prediction in low Pt bin 
(up to Pt < 35), but prefer the high order predictions (dyrad, madgraph, powheg, fewz, and PYTHiA+ljet 
) in Pt > 35. At low Pt, the measured values of A^ and A^ are well discribed by the qq — )■ ^* jZg annihilation 

function Aq — A2 ~ pi^M'^ ■ high Pt, the larger values show that both the annihilation and Compton 
process contribute to the cross section. Our results are in agreement with the fixed order perturbation theory 
calculations including dyrad, madgraph, powheg, fewz, and PYTHiA+ljet. In addition, we confirm the 
Lam- Tung relation (Aq ~ A2) which is only valid for vector gluon theory. The average of Aq — A2 is measured 
to be 0.017 ± 0.023, which is consistent with zero within the uncertainty. 

The measured A3 and A4 parameters are relatively flat in Z Pt. The central value of A3 parameter is 
slightly going down in Pt, but consistent with zero within one deviation of the uncertainty. The A4 parameter 
determines the forward and backward asymmetry of Drell-Yan process which has a mass dependence. Here, we 
measure A4 integrating of the mass range, 66 < Mee < llGGeV/c^. The Pr dependence of the average value 
of A4 for the mass range, 66 < Mee < 116 GeV/c^, is in agreement with the predictions of all the models. The 
A4 parameter is sensitive to Weinberg angle, sindyy. To compare the measured A4 with the predictions, we 
generate all predictions setting sin 9^ — 0.232. 




Figure 4: The measured Aq and A2 and predictions in Pt for the mass window, 66 < Mee < 116 GeV/c^. Aq and A2 
have a strong Pt dependence, the Lam- Tung relation {Ao ~ A2) which is only valid for vector gluon theory is confirmed. 
(The average of Aq — A2 is measured to be 0.017 ± 0.023.) For the measured Aq and A2, the uncertainties correspond 
to the total uncertainties {stat. ® syst.) 



6. Afh Measurement 

Afb is measured using the event counting method in mass. 

_ Ns^g{cose > 0) - NsigjcosO < 0) 

" Ns^3 (cos 9>0)+ Nsig (cos e<0) ^ ' 

where Nsig is the number of the signal events in the mass bin. The measured Aft using 4.1 fb~^ data is shown in 
Figure [7] and it is compared with pythia prediction. This measured A fb is unfolded using the response matrix 
(Rij) inversion method to get Ajb in the physics level. (true value) = P~^(i^(observation) - /3 (background)) 
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Figure 5: The measured ^3 and A4, and predictions in Pt for the mass window, 66 < Mee < 116 GeV/c^. A3 and A4, are 
relatively flat in Pt- For the measured A3 and At, the uncertainties correspond to the total uncertainties {stat. © syst.) 



The response matrix includes the detector smearing effect, the acceptance and efficiency {A x e) effect, and the 
final state photon radiation (FSR) effect. The response matrix is shown in Figure |6l 




Figure 6: The response matrix for unfolding Aft measurement. The response matrix includes the detector smearing 
effect. Axe, and FSR effect. 



The measured Afb is unfolded using the response matrix inversion and the unfolded Afb is shown in Figure 
[T] The unfolded Afb is compared to PYTHIA prediction and it has a good agreement with the prediction. 



7. Weinberg angle (sin^ 6w) vs. 



One of the angular coefficients, A4, is sensitive to the weak mixing angle, sin^ 9w- We translate A4 measure- 
ment into sin^ 9w using the various predictions to extract the weak mixing angle in 66 < Mee < 116 GeV/c^. 
The measured A4 in Pt is integrated over Pt and 66 < Mee < 116 GeV/c^. For the predictions, pythia, 
RESBOS, VBP, POWHEG, and FEWZ are used to extract sin^ 6w Figure [8] shows A4 vs. sin^ from various 
theory predictions. The theory band in Figure [5] includes the order diffdrence of the calculation and also PDF 
difference (CTEQ vs. MSTW). The variation in the predictions is assigned as the uncertainty of QCD model. 
The extracted sin^ 6^ = 0.2329 ± 0.0008 (A4 error);°;™J° (QCD). 
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Figure 7: The measured Afb in dielectron mass before and after unfolding using 4.1 fh~ data. The left plot shows the 
measured Aft, before unfolding and the right plot shows Afb after unfolding. For Afb before unfolding, the statistical 
uncertainty is only considered. For Afb after unfolding, the blue bars correspond to the statistical error and the red bars 
to the total error (stat. © syst.). The Afb measurement is compared to pythia prediction. 




8. Conclusion 

We present the results of the first measurement of the angular coefficients in the production of 7*/^ bosons at 
large transverse momenta, and the first test of the Lam- Tung relation at high transverse momentum using 2.1 
fh~^ data. We find a good agreement with the predictions of QCD fixed order perturbation theory, and with 
the Lam- Tung relation Aq ~ A2 . From measurement, we extract Weinberg angle, s\t?9w, using various 
theory predictions. We also present Afi, measurement in mass using 4.1 fh~^ data and the measured A^f, is 
compared to pythia prediction. The measured Afb has a good agreement with PYTHIA prediction. We expect 
that a comparison of these result with the measurement at LHC would provide additional tests of production 
mechanisms because the contribution of Compton process at LHC is expected to be larger than Tevatron. The 
angular measurement is published in Ref. [25|. 
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